Diffusive Domain Coarsening : Early Time Dynamics and Finite-Size Effects 
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We study diffusive dynamics of phase separation in a binary mixture, following critical quench, 
both in spatial dimensions d — 2 and d — 3. Particular focus in this work is to obtain information 
about effects of system size and correction to the growth law via appropriate application of finite-size 
scaling method to the results obtained from Kawasaki exchange Monte Carlo simulation of Ising 
model. Observations of only weak size effects and very small correction to scaling in the growth law 
are significant. The methods used in this work and information thus gathered will be of paramount 
importance in the study of kinetics of phase separation in fluids and other problems of growing 
length scale. We also provide detailed discussion on standard methods of understanding simulation 
results which may lead to inappropriate conclusions. 

PACS numbers: 64.60.Ht, 64.70.Ja 



I. INTRODUCTION 

When a homogeneous binary mixture (A+B) is 
quenched inside the miscibility gap, the system falls out of 
equilibrium and moves towards its new equilibrium state 
via formation and growth of domains rich in A- or B- par- 
ticles [lM3|. This coarsening of domains is a scaling phe- 
nomenon, e.g., two-point equal-time correlation function 
C(r, t), the structure factor S(k,t) and the domain size 
distribution function P(£d,t) obey the scaling relations 



C(r,t) = C(r/e(t)), 
S(k,t) = £ d S(k£(t)), 
P(£ d ,t) = l(t)-ip[l d /l(t)\, 



(1) 
(2) 
(3) 



where the average domain size £(t) increases with time (t) 
in a power-law fashion 



£(t) ~ r, 



(4) 



and C(x), S(y) and P(z) are scaling functions indepen- 
dent of £{t). In Eq. (j4j), the growth exponent a depends 
upon the transport mechanism. 

For diffusive growth, associating the rate of increase 
of £{t) with the chemical potential (/i) gradient, one can 
write [l| 



d£{t) 
dt 



IV" 



Ml 



£(t) 2 



(5) 



(7 being the A-B interfacial tension. Solution of Eq. J5j) 
gives a = 1/3, known as Lifshitz-Sloyozov (LS) law [4[. 
The LS behavior is the only growth law expected for 
phase-separating solid mixtures. However, for fluids and 
polymers, one expects faster growth at large length scales 
where hydrodynamic effects are dominant. For the latter, 
in d = 3, convective transport yields additional growth 
regimes Q with 



a = 1, £(t) <C £ in 
= 2/3, £{t) > £ in 



In Eq. (j6]), the inertial length £i n [— rj/ (per), r] and p being 
the shear viscosity and mass density] marks the crossover 
from a low- Reynold-number viscous hydrodynamic regime 
to an inertial regime. 

Diffusive domain coarsening in solid binary mixtures 
has been extensively studied via Ising model 



H 



j x SiSj] Si 

<ij> 



±1, J > 0, 



(7) 



prototype for a large class of critical phase transitions. 
Here one can identify the spin Si = +1(— 1) at lattice 
site i with an A-particle (B-particle). Note that < ij > 
in Eq. ([7]) stands for summation over only the nearest 
neighbors. One can also study the kinetics of phase sep- 
aration via dynamical equations which can be obtained 
from Ising models in mean field approximation by using a 
master equation approach [7|, [8| with Kawasaki exchange 
kinetics [9|. Upon coarse-graining, such equations lead to 
the Cahn-Hilliard (CH) equation 



dijj(r, t) 
dt 



V 2 m^t) + vW,t) - ^{r.t)], (8) 



(6) 



where ^(r, t) is a coarse-grained time-dependent local or- 
der parameter. Note that such continuum description 
could also be obtained in a phenomenological manner 
p], [loj using a coarse-grained Ginzburg-Landau(GL) free 
energy functional with the requirement of conservation of 
material. The CH equation with an added thermal noise 
is expected to be equivalent to Monte Carlo (MC) simu- 
lations [ll|, [l2| of kinetic Ising models. 

In Eq. (J5]), typical distance over which the order pa- 
rameter is coarse-grained is of the size of equilibrium cor- 
relation length, £. In situations when one is interested in 
studying the kinetics in the close vicinity of the critical 
point, without focusing on the dynamics at microscopic 
level, Eq. (|8j) is computationally very useful in achieving 
the asymptotic [13| behavior. However, for deep quenches 
one needs to incorporate higher order terms than are usu- 
ally used in the GL hamiltonian. Also at very low tem- 
perature, where £ is of the order of a lattice constant, CH 
equation would not provide information of a large effec- 
tive system size compared to the atomistic Ising model. 
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Particular focus of this work is to learn finite-size effects 
and dynamics at the early stage both of which have re- 
ceived much less attention as opposed to the identification 
of domain growth law in long-time limit, despite their ob- 
vious importance both fundamentally as well as techno- 
logically, e.g, in nano-science and technology. In view of 
that, we choose to revisit kinetics of phase separation in 
Ising model via MC simulations. 

While MC simulations have been used immensely in 
the understanding of non-equilibrium domain growth phe- 
nomena both with conserved |14l-l2l| and non-conserved 
[2ll-l23| order parameter, earlier studies of phase ordering 
in conserved Ising model with critical (50:50) composition 
reported [24], [H| estimates of a e [0.17,0.25], deviating 
drastically from the expected LS law. Even arguments in 
favor of logarithmic growth were proposed [26| . Note that 
these earlier reports were based on MC simulations for 
very short period of time where contamination of domain 
structures due to thermal noise might not have been taken 
care of, which could act as a source of significant error in 
the measurement of average domain size. 

Later, the discrepancy of the previous results with the 
expected LS behavior was understood to be due to strong 
corrections to scaling at early time. To account for this 
[l5| higher order terms in Eq. (j5j) were incorporated to 
write 



d£(t) _ d 
dt ~ 



C 2 



(9) 



which in the long time limit gives a solution oc t 1 / 3 , how- 
ever, would give rise to leading order correction linear in 
l/£(t) to the instantaneous exponent. Thus LS behavior 
will be observed only in the limit £(t — >• 00) — »> 00. In- 
deed, consistency with a linear correction was observed for 
50:50 binary mixture |15|, [l6| as well as for multicompo- 
nent mixtures [19]. Present work, however, convincingly 
demonstrates that the observation of LS value of the ex- 
ponent only in the asymptotic limit was misleading and 
presence of a time independent bare length in l(t) was 
responsible for the numerical data exhibiting such trend. 

Most of the studies till date, stressed on using large 
systems, with the anticipation of strong finite-size effects 
|27j combined with the expectation that the LS law will 
be realized only in the large £(t) limit. This strategy, of 
course, will prove to be useful when there is dynamical 
crossover as in domain coarsening in fluids [cf. Eq. (0] 
where the system size should be significantly larger than 
smallest characteristic length scale in a particular regime. 
However, consideration of arbitrarily large system sizes 
restricts the access of large time scale, particularly for 
molecular dynamics simulation of fluid phase separation 
[28l-[30|. It is worth mentioning that the typical system 
sizes authors consider these days contain number of lattice 
sites or particles of the order of million, which is too large, 
even for present day computers, to access long time scale 
that often is a necessity. Such choice of large systems, in 
addition to the anticipation of strong finite-size effect, was 
often motivated by the expectation of better self- averaging 
[27( | . However, our experience contradicts the latter and 



is rather consistent with the previous works [22] reporting 
lack of self averaging. Thus a judicial choice of system sizes 
is very crucial for such problems [31] which in turn requires 
appropriate knowledge of finite-size effects [32| . While re- 
cent focus has been in more complicated systems [33l-[40| 
with realistic interactions and physical boundary condi- 
tions, many such basic information as discussed above are 
lacking even in a situation as simple as Ising systems. 

In this work we have undertaken a comprehensive study 
to learn about the finite-size effects in domain coarsening 
in Ising model with conserved order parameter dynamics 
and understand the behavior of growth exponent as a func- 
tion of time, via application of finite-size scaling method 
[4lL |42|, both in space dimensions d = 2 and d = 3. De- 
spite its glorious popularity in equilibrium critical phe- 
nomena, finite- size scaling method has been only rarely 
[20I , l2ll K32j used in non-equilibrium processes. In this pa- 
per, we exploit this method appropriately in the context 
of diffusive phase separation kinetics to show that for crit- 
ical quench the LS value of a sets in very early and the 
effect of size is very small. 

This paper is organized as follows. In Sec Jill we 
describe the details of simulation and finite-size scaling 
method. Results for both d = 2 and d = 3 are presented 
in Sec. IHIi while Sec. IIVI summarizes the paper with a 
discussion of future possibilities in this direction. 



II. METHODS 

A. Details of Simulation and Calculation of 
observables 

In the MC simulation of Ising model, the conserved 
order-parameter dynamics, where composition of up (A 
particle) and down (B particle) spins remains fixed dur- 
ing the entire evolution, is implemented via standard 
Kawasaki exchange mechanism [9] where interchange of 
positions between a randomly chosen pair of nearest neigh- 
bor (nn) spins consists a trial move. A move is accepted or 
rejected according to standard Metropolis algorithm [TT| . 
One MC step (MCS) consists of exchange trials over L d 
pairs of spins. Periodic boundary conditions were applied 
in all directions. 



FIG. 1: Left panel: Snapshot of a 2 — d Ising model at T — 
0.85T C obtained from the Monte Carlo simulation via Kawasaki 
exchange kinetics, for L = 64 at t = 5 x 10 3 MCS. Right 
panel: Same snapshot after removing the noise via the exercise 
described in the text. A-particles are marked by black dots 
whereas B-particles are unmarked. 
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Note that with the increase of temperature, accurate 
measurement of average domain size becomes difficult due 
to the presence of noisy clusters of the size of £(T). On the 
other hand, at very low temperature growth is hampered 
by metastability. To avoid the latter problem, we have set 
the temperature towards the higher side and calculated all 
the physical quantities from pure domain morphology af- 
ter eliminating the thermal noise via a majority spin rule. 
There a spin at a lattice site i was replaced by the sign of 
the majority of the spins sitting at i and nn of i (depending 
upon the noise level i.e., average size of noise clusters, ex- 
tension to distant neighbors may also become necessary). 
In Fig. [T]we demonstrate such filtering process for a rather 
high temperature. The left panel corresponds to the orig- 
inal snapshot from the MC simulation on 2 — d square 
lattice at T = 0.85T C with L = 64 at t = 5 x 10 3 MCS. 
One can appreciate that presence of substantial noise el- 
ements here would give rise to smaller value of £(t) than 
the actual. The right panel of the figure shows the snap- 
shot with pure domain morphology obtained after imple- 
menting the noise removal exercise described above. Of 
course, one should be careful that too many such itera- 
tions or consideration of very distant neighbors may alter 
the basic structure. However, in the present case, no such 
deformation has taken place. All the quantities in our 
simulation were calculated by using snapshots with such 
pure domain structure. 

In Fig. [2fa) we present the scaling plots of domain size 
distribution function, viz., plots of £{t)P(£d,t) vs £d/£(t) 
where £(t) was calculated from the first moment of the 
normalized distribution P(£d,t) as 



£{t) 



l d £ d P(£ d ,t), 



(10) 



with length £& being the separation between two successive 
interfaces in a>, y- or z- directions. Figs. EJb) and [2fc) 
show the scaling plots of correlation function C(r,t) and 
its Fourier transform S(k,t), in accordance with Eq. (pQ) 
and (2), where C(r, t) was calculated as 



C(r,t) =<SiSj >; 



r = \i 



(ii) 



Note that these scaling plots for all the quantities were 
obtained by using the values of £{t) obtained from Eq. 
([TO]) . Of course, independently £(t) could be calculated 
from the decay of C(r, t) as well the first moment of S(k, t) 
as 
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FIG. 2: Scaling plots of (a) domain size distribution P(£d,t), 
(b) correlation function C(r, t) and (c) structure factor S(k,t), 
from different times, as indicated, for the system in Fig. [T] The 
data were averaged over 50 independent initial configurations. 
Note that in all the cases, £{t) used, was calculated using Eq. 
(|10p . While the solid line in (a) verifies exponential decay of 
the tail, the one in (c) corresponds to the Porod tail. 



C(r = l(t),t) =ft, 



(12) 



and 



£(t) = 



Jdk kS(k)' 



(13) 



When calculated from a completely noise- free morphol- 
ogy, all the above mentioned methods for the calculation of 
£{t) must give results proportional to each other. When h 
is set to a rather small value, particularly when the decay 
length is larger than the average size of the noisy clusters, 



calculation of £(t) from Eq. ([T2]) is not expected to be 
affected much by the presence of noise. The same applies 
for Eq. (fT3|) . However, when calculated via Eq. ([TO]) , 
either the distribution up-to the length of average noise 
size should be appropriately modified or noise should be 
completely eliminated. The latter strategy is more appro- 
priate since it gives better shape to the all form functions. 
In our calculation, in Eq. ([T2]h h will correspond to first 
zero of C(r, t). 
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All the results presented in Fig. [2] are obtained from 
pure domain morphology and the nice data collapse ob- 
tained in each case using the measure of l(t) from Eq. 
([TO]) speaks for the equivalence of all the definitions, Eqs. 
(fTU|) , (fT2j) and (fT3|) . The linear behavior of the tail region 
in (a) on a semi-log plot is consistent with an exponential 
decay of P(£d,t). Here the noisy look (oscillatory behav- 
ior) at late time or large domain size limit (which was 
also observed in other recent studies [43|, is due to 
lack of statistics when l(t) becomes of the order of the 
system size. We refer to it as a secondary finite-size effect 
which does not affect the growth dynamics severely. On 
the other hand, the linear look of large wave vector (k) 
data in (c) confirms the Porod law 



#(M)~£;- (d+n) . 



(14) 



Note that in the present case d = 2 and n = 1 (symmetry 
of the order parameter). It is worth mentioning that one 
would have observed a slower decay of the structure factor 
had the noise not been removed. 
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FIG. 3: Average domain size is plotted on a log-scale as a func- 
tion of time t. Different symbols correspond to calculation of 
£(t) from different quantities : circles from P(£d,t), squares 
from first zero-crossing of C(r, £), diamonds from the first mo- 
ment of S(k,t). Results presented were obtained from pure 
domain structure as demonstrated in the right panel of Fig. [1] 
with L 2 — 128 2 and T = 0.85T C . The solid line corresponds to 
the theoretically expected t 1//3 behavior. 



In Fig. [3] we present the length scale results obtained 
from all the above mentioned methods on log scale, all of 
which look proportional to each other as was also clear 
from the exercise of Fig. [2j The data at late times look 
consistent with the expected exponent a — 1/3. Note 
that if the temperature is sufficiently close to T c , for long 
enough time the noise might not have equilibrated to the 
equilibrium value inside the true domains. In such a situ- 
ation, there are two growing length scales in the problem 
which may give rise to misleading value of the exponent 
if the noise is not eliminated and range of fitting is small. 
Indeed a fitting of the data, obtained from original snap- 
shots (not shown) to the form 



in the range [0,20000] MCS gives a = [0.15,0.25] (the 
value being larger when £(t) is calculated from ([12]) or 
([T3]) ) which is consistent with earlier reports [24], |25j. On 
the other hand a similar fitting to the data obtained after 
removing the noise gives a = [0.3, 0.34] and within statis- 
tical deviation, does not depend upon the range of fitting. 
This latter result is already suggestive of absence of strong 
correction to the scaling. However, since data fitting is al- 
ways not a very reliable exercise as will be discussed later, 
to further substantiate the claim about small correction 
to scaling, we take the route of finite-size scaling analysis 
that will also be useful in quantifying the finite-size effect. 



B. Formulation of Finite-Size Scaling 

In equilibrium critical phenomena, the singularity of 
a quantity Z is characterized in terms of e = \T — T c |, 
temperature deviation from the critical point, as 



Z « Z e z « Z'& 



-z I V . 



z' 



y t z l v 

z o4o 



where the correlation length £ grows as 



(16) 



(17) 



with z and v being the critical exponents. However, for 
finite values of L any critical enhancement is restricted 
and Z is smooth and analytic. Such finite-size effects may 
appear as a difficulty in understanding results from com- 
puter simulations. However, this problem can be tackled 
by writing down finite-size scaling ansatz, thus accounting 
for the size effect, as 



Z * Y(x)e z = Y\x)C 



z I V 



(18) 



In Eq. (fT8|) . Y(x) is the finite-size scaling function that 
depends upon the scaled variable x = L/£ and is indepen- 
dent of system size. Note that Y and Y' differ by a factor 
originating from different amplitudes Zq and Z' Q used in 
Eq. ([T6]) . In further discussion, however, we will remove 
primes from both Zq and Y and distinction can be derived 
from whether the scaling forms are written in terms of e 
or £. 

At this stage, it is important to ask about the behavior 
of Y as a function of x. While for static quantities, such 
question is already addressed, for dynamics, where the 
finite-size effects are found [45| to be much stronger, there 
is no appropriate understanding of the variation of Y{x). 
Nevertheless, one can write down the following limiting 
behaviors: 



for x ->> (e ->> 0; L < oo), Y(x) ~ x 
such that Z is finite at criticality : 
Z - L~ z/v . 



-zjv 



(19) 



(20) 



£(t) = C + At a , 



(15) 



Eq. ([20]) , when compared with Eq. ([T6]h is consistent 
with the fact that at criticality £ is the only important 
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length in the problem and L must scale with £. Keeping 
this important fact in mind, in fact, one can write ([T8]) as 



Y{x)L 



-z jv 



On the other hand, 



for x oo (L — >• oo, e > 0), Y"(x) = Zq, 



(21) 



(22) 



so that Eq. ([T6]) is recovered in the thermodynamic limit. 

With the knowledge of Eq. ([20]) can be used to esti- 
mate z by calculating Z at T c for various system sizes. A 
better strategy however is to study Z at finite-size critical 
points T C L , such that 



-z/v . 



-l/i/ 



(23) 



though true meaning of a critical point can be assigned to 
T C L only in the limit L — >> oo. 

This general discussion about finite-size scaling method 
could be used to construct the similar formalism for non- 
equilibrium domain coarsening problem by identifying l(t) 
with £ and 1/t with e. In the present problem, L should 
scale with £{t), more precisely £ max ~ £(£), where f max (L) 
is the equilibrium domain size for a system of size L. In 
Fig. |4ja) we show plots of £(t) vs £ for various different 
values of L, in d = 2. The flat regions in the plots at late 
times correspond to ^ max . 

At this stage, we would like to quantify the domain 
growth as 



£(t>) 



At" 



(24) 



where to is the average cluster size when the system be- 
comes unstable to fluctuations at time to since the quench. 
Of course our measurement of time starts from there, i.e. 
t' = t — to. Note that we do not assign a meaning of 
domain size to this quantity and this should be treated 
in a manner similar to a background quantity in critical 
phenomena that appears from small length fluctuations 
whose temperature variation is usually neglected. Hav- 
ing said that, scaling part in Eq. (|2^|) is only At /a . Of 
course, when t(t') is significantly large, subtraction of the 
microscopic length to does not bring in noticeable differ- 
ence. However, in computer simulations, where one deals 
with small systems, presence of to can bring in completely 
different conclusions. 

Eq. is valid only in absence of any finite-size effect. 
For a finite 4ax(^), analogous to ([H]), one can write down 
the scaling ansatz as 



where now 



£{t') -£o = Y(x)t' 1 



(25) 



(26) 



is the scaling variable. Both in Eqs. ([25]) and ([26]) . to is 
subtracted to deal with the scaling parts only. By observ- 
ing (USD, CEO) and fl22j) as well as (JM]), (gSJ) and ([26]), one 
can arrive at the limiting forms of Y(x) as 

Y{x) w x, for x -> 0(t f -> oo,^ max < oo) (27) 
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FIG. 4: (a) Plot of average domain size £(£), obtained from 
Eq. ([TQ]h for the 2 — d Ising model, for different system sizes 
(indicated on the figure) at T — 0.6T C . Definitions of ^ ma x 
and tf^ are demonstrated. Data for L — 32 and L — 64 were 
averaged over 1000 independent initial configurations whereas 
only 40 different initial realizations were used for L — 128. 
Note that all subsequent results in this paper are obtained at 
the same temperature as this, (b) Demonstration of the scaling 
behavior (|29[) in d — 2 and 3. 



and 



Y(x) = A, for x — > oo(t' < oo, 



oo). (28) 



Of course, it would again be interesting to learn about the 
full form of Y(x). 

Also, analogous to T C L in critical phenomena, one can 
define a finite size equilibration time t^ q that is needed to 
reach ^ max , as demonstrated in Fig[4]^a). Then one can 
write down a scaling equation analogous to ([23]) as 



[4 



1/3 



(29) 



This scaling behavior is demonstrated in Fig. H(b) where 
we show plots of l max vs t^ q on a log scale, in both d = 2 
and d = 3. Consistency of the simulation data with the 
solid line of the form ([29]) confirms the validity of this 
approach. Note that in this figure, we did not subtract £q 
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and the corresponding microscopic time from the abscissa. 
As will be seen later, £ max for the systems considered here 
are significantly larger than so one does not expect a 
big difference after subtracting. Eq. ([29]) , of course, is a 
statement of the fact that dynamic critical exponent for 
this model is 3. 



III. SIMULATION RESULTS AND ANALYSIS 

Having set the methodology in place, in this section 
we present results from MC simulation of Kawasaki-Ising 
model in d — 2 and 3, combined with the finite-size scaling 
analysis. 



A. Results in d=2 

In Fig. [5] we present snapshots during the evolution 
of an Ising system, starting from a 50 : 50 random mix- 
ture of up and down spins, obtained via MC simulation 
at temperature T = 0.6T C . The times at which the shots 
were taken are mentioned on the figure. While the last 
snapshot corresponds to a situation when A and B phases 
are completely separated, the one at t = 3.5 x 10 5 MCS 
represents the situation when finite-size effect began to 
enter, which will be clear from subsequent discussion. 





3.5X10 mcs 



10 MCS 





FIG. 5: Evolution snapshots from different times, as indicated, 
for the Kawasaki-Ising model in d — 2 at T — 0.6T C . The last 
snapshot corresponds to a completely equilibrated configura- 
tion. 

From Fig. [4(a), it is already evident that finite-size ef- 
fect is rather small. However, for a quantitative statement 
and to gain detail information about the growth expo- 
nent, more sophisticated analysis is called for. Following 
the discussion in the previous section, in Fig. [6] we plot 
Y = [£{t') - £ }/t' a as a function of x/(x + x ). Note that 



xo was introduced to see behavior of Y properly both for 
small and large x. For convenience we set it to 5. In this 
exercise we have varied a and £q (or the microscopic time 
to associated with this length) to get optimum collapse of 
data from different system sizes. In Fig. [5(a), where £{t) 
is being used from Eq. ([TO]) , the optimum data collapse 
is obtained for £q ~ 4a (average cluster size after 20 MCS 
since quench), a being the lattice spacing and a ~ 0.33. 
Similar exercise when £{t) is being obtained from Eq. ([T2]h 
as shown in Fig. 15(b), gives a ~ 0.35 and £q there corre- 
sponds to the same number of MCS after quench. Note 
that £o in our analysis is a bare length independent of 
time and the scaling behavior ([25]) will be obtained when 
this is chosen appropriately. These numbers, as expected, 
provided a constant value of Y(x) in the region unaffected 
due to finite system size, which should be identified with 
the growth amplitude A for which we quote 0.29 ±0.01 (cf. 
Fig. [5(a)). The arrows in Fig. [5] marks the location where 
Y(x) starts deviating from its constant value. The sharp 
nature of the crossover is indicative of only small size effect 
which we quantify from the location of the arrow marks 
as 



£(t) = (0.75 ±0.05)4 



(30) 



Of course, this value is significantly large compared to ear- 
lier understanding and expectation. Note that the snap- 
shot at t = 3.5 x 10 5 MCS in Fig. [5] corresponds to this 
length. 

In an attempt to learn the full form of Y(x), we con- 
struct the following functional form 



Y(x) 



Ax 



x ± l/(p ± qx@) ' 



(31) 



that has limiting behaviors consistent with ([27|) and ([28]) . 
The continuous lines in Figs. [5(a) and (b) are fits to the 
form (I3T1) with 



A ~ 0.29, q ~ 6400, j3 = 4 



(32) 



and 



A ~ 0.14, p~7, q~ 13700, /3 = 4, (33) 
thus have the convergence 



(x -> oo) Y{x) « A[l - fx~ n ];n = 5. 



(34) 



Of course, possibility of an exponential correction cannot 
be ruled out. This may be compared with much slower 
convergence of such function in dynamic critical phenom- 
ena [45[. Note that the understanding of finite-size effect 
in both equilibrium and non-equilibrium dynamics is a 
non-trivial task and significant attention is called for. 

To take a direct view of what happens after the cor- 
rective measure, in terms of subtraction of £o, has been 
taken, in Fig. [3 we plot [£{t') - £ }~ 3 vs l/t' and £(t)~ 3 
vs 1/t for L = 64. A log-scale was used to bring visibil- 
ity to a wide range of data. The linear behavior of the 
data after subtracting £q : starting from very early time 
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0.1 



t — i — i — i — | — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 

i 6b 6 to a : 




\\£(t)=0.75£ max 



o L = 32 
□ 64 
o 128 
— Eq.(31) 



i ■ ■ ■ ■ i ■ ■ ■ ■ i ■ ■ ■ ■ i 



(a) 
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0.04 
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£(t) = 0.75W ° L = = f A 

o = 128 
— Eq. (31) 
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x/(x + x ) 



0.5 



FIG. 6: (a) Finite-size scaling plot of Y, with £o = 3.6 lattice 
constants (after 20 MCS from the quench time) and a ~ 0.33, 
as a function of x/(x + xo); xo = 5. The continuous curve is a 
fit to Eq. (|31[) with the best fit parameters mentioned in the 
text. The arrow roughly marks the appearance of finite-size 
effect. Note that £{t) data used here, were obtained from Eq. 
([TO]) , (b) Same as (a) but £{t) obtained from the first zero- 
crossing of C(r,i) [cf. Eq. (p"2]) ]. In this case £o — 2.7 lattice 
constants (at 20 MCS from quench) and a c± 0.35. 
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FIG. 7: Plot of [£(t') - £o}~ 3 vs 1/t' and £(t)~ 3 vs 1/t for 
L 2 = 64 2 , with £(t) being calculated from Eq. (TJJ. The 
continuous line has slope 39 = 1/A 3 . 
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FIG. 8: Scaling plot of C(r,t) at T = 0.6T C . Note that £{t) 
was obtained using Eq. ([TO]) . 



justifies the introduction of £q again. The continuous line 
there is a plot of the form Ax with A ~ 39 = 1/A 3 . On 
the other hand, notice the strong curvature when £$ is not 
subtracted. The dashed lines marked by 1/t and ( j) 3 ^ 5 on 
this figure corresponds to £(t) ~ t 1 ! 3 and t 1 ^ respectively. 
Thus, when £q is not appropriately subtracted, only ob- 
serving the trend on a log-log plot one may be misled to 
conclude that there is gradual crossover from one regime 
to the other. See Ref [46| for a discussion on earlier belief 
about a crossover from t 1 / 4 to t 1 ^ 3 . Note that the exercise 
here as well as one in Fig. [6j where Y is very fiat from 
very early time all the way to the moment when finite- 
size effect enters, are already indicative of absence of any 
strong corrections to scaling. 

Before moving ahead for another elegant proof of the 
evidence for the absence of negligible corrections to scal- 
ing, we pass by showing the scaling plot of C(r, t) in Fig. [8] 
where good quality data collapse is obtained starting from 
very beginning till t = 3 x 10 5 MCS when the finite-size 
effect begins. Next we introduce a length £ s to write 



£'(t') = £(t')-£ s = [£o-e s } + At' a , 
and calculate the instantaneous exponent [15 
_ rf[lnf (f)] 



d[lnt'} 



to obtain 



1 



£>(t>) 



(35) 



(36) 



(37) 



According to Eq. ([37]) . when oti is plotted as a function 
of l/£'(t'\ for £'{t') > 0, one expects linear behavior with 
a ^-intercept equal to a. Fig. [9fa) shows such plots for 
£ s = 0.0,3.6, and 5.0, as indicated. The dashed lines have 
^/-intercept a = 1/3 and slopes 



to -i* 



(38) 



The consistency of the actual data with the dashed lines 
is remarkable, particularly the behavior of for £ s = 3.6, 
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again speaks for the choice of £o and strongly indicates 
that the LS scaling regime is realized very early. In Fig. 
E^b) we present results with £ s = 3.6 for various system 
sizes 1? = 16 2 ,32 2 and 64 2 . In all the cases, oscil- 
lates around 1/3. This observation, using a system size as 
small as L 2 = 16 2 , stresses against unnecessary attempt 
to simulate larger systems. 



0.45 



OLA 




0.36 

0.35 

OLi 

0.34 

0.33 

0.32 

0.31 



~i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 



L=64 L = 32 - % -£=16" 
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i/i'{t') 



FIG. 9: (a) Plot of instantaneous exponent ol% as a function 
of l/i'{t') for three different choices of i s as indicated, with 
L 2 = 64 2 . The dashed straight lines have slopes —1.19,0 and 
0.49, respectively, (b) Plot of on vs l/£'(t') for i s = 3.6 and 
L 2 = 16 2 ,32 2 and 64 2 . In both (a) and (b) the arrows on the 
ordinate mark the value a = 1/3. Note that £{t) was calculated 
from Eq. (HOT) . 



of negligible correction to the exponent, starting from the 
very early time, is consistent with the growing evidence 
[HI, l49l th at Tolman length is absent in a symmetrical 
model [50j where the leading correction is of higher order. 
Also, small corrections that may be present, coming from 
the curvature dependence of the kinetic pre-factor in Eq. 
([5]) is beyond the accuracy of data in the present work. 
On the other hand, for 50 : 50 composition, since the 
domain boundaries are essentially flat starting from very 
early time, any curvature dependence is expected to be 
absent. Thus, we conclude that this misunderstanding 
about the strong time dependence in a was due to the 
presence of a time independent length £0 in £(t) which our 
analysis subtracts out in appropriate way. 



5 X 10 4 MCS 
*** 




MCS 



FIG. 



10: Evolution snapshots from different times for 3 — d 
Ising model with L 3 = 32 3 and T = 0.6T C . A and B particles 
are marked black and grey respectively. 



This result is in strong disagreement with the earlier 
[IH understanding of domain coarsening in 2 — d conserved 
Ising model for critical quench that a is strongly time de- 
pendent and the LS value is recovered only asymptotically 
as t(t — > 00) — ^ 00. The route to this finite-time correc- 
tion was thought to be an additional term oc l/£(t) 3 in 
Eq. ([5|) [cf. Eq. ©], accounting for an enhanced inter- 
face conductivity. Note that a term ex l/£(t) 3 could also 
be motivated by introducing a curvature dependence in a 



a[£(t)} 



<t(oo) 

1 7 2S ■ 



(39) 



S being the Tolman length [47]. However, our observation 



B. Results in d=3 

In this subsection we turn our attention to the kinetics 
of phase separation in d = 3. Fig. [T0l shows 3 — d snap- 
shots of the time evolution of Kawaski-Ising model at four 
different times as indicated on the figure where the last 
snapshot is clearly seen to have been equilibrated. Anal- 
ogous to d = 2, all results presented here were obtained 
at T = 0.6T C , with T c = A.51k B T/J in this case, and the 
composition was chosen to be 50 : 50 as well. 

In Fig. [TTTa) we present direct plots of £(t) as a function 
of t, for L 3 = 16 3 ,32 3 and 64 3 where £(t) was calculated 
from Eq. ([T2]) . Again, the finite-size effects look to be 
small. In Fig. Qjjb) we present a plot of Y(x), using the 
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; £{t) = 0.7€ max 




o L =32 
□ 64 
— Eq. (31) 



(b): 
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FIG. 11: (a) Plot of £(t), obtained from the first zero crossing 
of C(r,t), vs t, for the systems L 3 = 16 3 ,32 3 and 64 3 . (b) 
Finite-size scaling plot of Y, for the data presented in (a), vs 
x/(x-\-xo) with xo = 5. Here to = 10 (10 MCS from the quench) 
and a = 0.35. Appearance of finite size effect, obtained from 
the arrow mark, is estimated to be at £{t) ~ 0.7<£ ma x, in close 
agreement with the one for d — 2. 



data in Fig. ITTTa), as a function of x/(x + xo);xo = 5. 
Best data collapse in this case was obtained for £q — 2.5 
(10 MCS after the quench, note that corresponding value 
of £ from Eq. JTGJ) is 3.0 and a ~ 0.315) and a ~ 0.35. 
Very flat behavior of Y(x), starting from the beginning 
again speaks for absence of any strong correction to the 
growth law. However, compared to the d = 2 case, one 
may expect slightly stronger correction here because of 
the inherent curvature present in the cylinder like domain 
objects as opposed to the stripe like structures in d = 2. 
Possibly because of that we could not obtain good collapse 
of data from L 3 = 16 3 on top of the ones presented, since 
the whole data set for L 3 = 16 3 is from very early time 
and suffers from corrections. Onset of finite size effect, as 
obtained from the arrow mark where Y(x) deviates from 
the fiat behavior, is in quantitative agreement with the 2- 
dimensional situation, as quoted in Eq. ([30]) . Here also the 
third snapshot in Fig. [10] (at t = 10 5 MCS) is presented 
at this onset. A fitting, shown by the continuous line, to 



the form (HIT]) [A = 0.24, p ~ 4, q ~ 13050 andj3^= 4], 
again is consistent with asymptotic convergence 



OLA 




£ s = 2.5 - 



0.2 0.3 0.4 
l/*'(f) 



FIG. 12: Plot of instantaneous exponent on vs l/£'(t') with 
three different values of £ s = 0, 2.5 and 5. The dashed lines 
correspond to a — 0.34. The arrow on the ordinate marks the 
value a — 1/3. 

In Fig. [12j we present instantaneous exponent as a 
function of l/£'(t') for L 3 = 64 3 and three choices of £ s as 
indicated. In all the cases, the exponent fluctuates around 
mean value 0.34. Note that a estimated from S(k,t) and 
P{£d,t) are slightly higher and lower, respectively, com- 
pared to the one presented. 

The appearance of growing oscillation in c^, seen in 
Figs. [9] and [12j around the mean value was also pointed 
out by Shinozaki and Oono [27[. In a finite system, as 
time increases, for an extended period of time two large 
neighboring domains of same sign may not merge, thus 
lowering the value of a. After a long time when two large 
domains merge, brings in drastic enhancement. This char- 
acter is in fact visible in the direct plot of £{t) vs t at late 
times [cf. L = 128 in Fig. gfa) and L = 64 in Fig. flTTa)]. 
Note that this oscillation could be a route to an error if 
one obtains a from least square fitting without choosing 
the range appropriately. Finally, it will be interesting to 
know the temperature dependence of £q and amplitude A 
as well as of finite size effects. All these, however, we leave 
out for future work. 



IV. SUMMARY 

This paper contains comprehensive study of domain 
coarsening in a phase separating system with diffusive dy- 
namics in d — 2 and d — 3. Various different ways of 
analysis give results for growth law consistent with the 
expected LS exponent a = 1/3. As opposed to the ear- 
lier understanding, correction appears to be very weak, 
thus LS scaling behavior being realized very early. Small 
primary finite-size effect is a welcome message which is 
suggestive of avoiding large systems, rather focusing on 
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accessing long time scale which often is necessary for sys- 
tems exhibiting multiple scaling regimes. 

Our observation should be contrasted with an earlier 
study of Heermann, Yixue and Binder [20| that reports 
very strong finite-size effect. However, this latter study 
was based on an extremely off-critical composition and 
should not be considered to have general validity. Note 
that due to the expected presence of correction in such 
off-critical composition, where droplet like structures form 
with finite radius of curvature at early time, the analysis is 
more difficult. Also, one should be prepared to encounter 
stronger size effect in more complicated situations, e.g., 
systems exhibiting anisotropic patterns (3314371 l39l Eof . 

One may of course ask if the small finite-size effect 
observed for diffusive dynamics is also valid for kinetics of 
phase separation in fluids. A comprehensive study in that 
direction, for both binary and single-component fluids, is 
in progress. Such studies are important since brute force 
method of simulating very large systems, particularly for 
the study of fluid phase separation via MD simulation, is 
not often helpful to access long time scales even with the 
present day high speed computers and thus may not bring 
very conclusive understanding. 

A deeper understanding of £o requires further study, 
particularly, how the system is led to instability is a fun- 
damental question to be asked. Even though scaling cor- 
rections appear to be negligible for critical quench due to 



the fiat nature of the domain boundaries, one expects cor- 
rections for off-critical composition. This expected correc- 
tion coming from surface tension should be of higher order 
than linear for a symmetric model. On the other hand, 
it would be interesting to learn about the leading order 
correction coming from the kinetic pre-factor. 

Finally, we expect the observation, understanding and 
finite-size scaling technique used in this work to find rel- 
evance in other systems exhibiting growing length scales, 
e.g., ordering in ferromagnets, surface growth, clustering 
in cooling granular gas, dynamic heterogeneity in glasses, 
etc. In line of this work many earlier studies on domain 
coarsening may need to be revisited for better understand- 
ing which was not possible because of lack of reliable meth- 
ods of analysis. 
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